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We present the results of a search for short-duration gravitational-wave bursts associated with 39 
gamma-ray bursts (GRBs) detected by gamma-ray satellite experiments during LIGO's S2, S3, and 
S4 science runs. The search involves calculating the crosscorrelation between two interferometer 
data streams surrounding the GRB trigger time. We search for associated gravitational radiation 
from single GRBs, and also apply statistical tests to search for a gravitational-wave signature asso- 
ciated with the whole sample. For the sample examined, we find no evidence for the association of 
gravitational radiation with GRBs, either on a single-GRB basis or on a statistical basis. Simulating 
gravitational- wave bursts with sine-gaussian waveforms, we set upper limits on the root-sum-squarc 
of the gravitational-wave strain amplitude of such waveforms at the times of the GRB triggers. We 
also demonstrate how a sample of several GRBs can be used collectively to set constraints on popu- 
lation models. The small number of GRBs and the significant change in sensitivity of the detectors 
over the three runs, however, limits the usefulness of a population study for the S2, S3, and S4 runs. 
Finally, we discuss prospects for the search sensitivity for the ongoing S5 run, and beyond for the 
next generation of detectors. 



I. INTRODUCTION 

It has been over three decades since gamma-ray bursts 
(GRBs) were first detected by the Vela satellites Dur- 
ing the 1990s, when the Burst and Transient Source Ex- 
periment (BATSE) and BeppoSAX were m oper- 
ation, important discoveries and observations relating to 
GRBs were made, such as the isotropic distribution of 
GRBs 0; the bimodal distribution of burst durations, 
suggesting long and short classes of GRBs pjjdetections 
of the first x-ray [||, optical 0, and radio M, counter- 
parts; the first redshift measurements 0, [H E| j and the 
first hints of the association of long-duration GRBs with 
core-collapse supernovae [H, 0, 03 ■ Today, important 
questions about GRB progenitors, emission mechanisms 
and geometry linger, and observations made by the cur- 
rent generation of gamma-ray satellite experiments such 
as Swift [H, HETE-2 0, INTEGRAL and others 
continue to provide new and exciting information which 
help us answer these questions and better understand the 
origin and physics of these astrophysical objects. 

Currently favored models of GRB progenitors are core- 
collapse supernovae for long-duration GRBs [HI], and 
neutron star-neutron star (NS-NS) or neutron star-black 
hole (NS-BH) mergers for short-duration GRBs 
These models and the division into two classes of progen- 
itors are supported by observations of supernovae asso- 
ciated with long-duration GRBs pj, El Q HH, H3] and, 
more recently, observations of afterglows and identifica- 
tion of host galaxies for short-duration GRBs [23J, \2A\ |25j, 
|26| . The end result in either scenario is the formation 
of a stellar-mass black hole 27\ and, in either scenario, 



theory predicts the emission of gravitational radiation. 
In the former case, gravitational waves would result from 
the collapse of a massive star's core, while in the latter 
case, gravitational radiation would result from the inspi- 
ral, merger, and ringdown phases of the coalescence. Re- 
cently, there has been an observation-driven suggestion 
of a third class of GRBs which could include both short- 
and long-duration GRBs (28| , but more observations are 
needed to support this suggestion. 

Due to the expected evolution of the proposed progen- 
itors, the redshift distribution of long-duration GRBs is 
thought to follow the star formation rate of the Universe 
[H (3^1 , and recent redshift measurements tend to sup- 
port this model, with the measured GRB redshift dis- 
tribution peaking at z > 1 31]. Long-duration GRBs 
have also been associated exclusively with late-type star- 
forming host galaxies [HJ ■ On the other hand, the recent 
observations of x-ray and optical afterglows from a few 
short-duration bursts seem to suggest that these GRBs 
are located at lower redshifts relative to long-duration 
GRBs [H, HI] , and that short bursts are found in a mix- 
ture of galaxy types, including elliptical galaxies, which 
have older stellar populations. All of these observations 
are consistent with the currently favored models of GRB 
progenitors. Although a large fraction of GRBs are too 
distant for any associated GW signals to be detected by 
LIGO, it is plausible that a small fraction occur at closer 
distances. For example, a redshift of z = 0.0085, or a dis- 
tance of 35 Mpc, has been associated with long-duration 
burst/supernova GRB 980425/SN 1998bw [Hf. It is not 
unreasonable to expect that a few GRBs with no mea- 
sured redshifts could have been located relatively nearby 
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as well. For short-duration GRBs, the recent redshift ob- 
servations have led to fairly optimistic estimates [34], HH 
for an associated GW observation in an extended LIGO 
science run. 

In this paper, we present the results of a search for 
short-duration gravitational-wave bursts (GWBs) asso- 
ciated with 39 GRBs that were detected by gamma-ray 
satellite experiments on dates when the S2, S3, and S4 
science runs of the Laser Interferometer Gravitational- 
Wave Observatory (LIGO) were in progress. Although 
the theoretical shapes of the GW burst signals result- 
ing from the two progenitor scenarios are not known, 
many models predict that the GW signals would be 
of short duration, ranging from ~ 1 ms to ~ 100 ms 
[36L l37l [HI, [H, [40] • The search method presented here 
targets such short-duration signals, and calculates the 
crosscorrelation between two LIGO interferometer data 
streams to look for these signals. A crosscorrelation- 
bascd method efficiently suppresses uncorrelated tran- 
sient noise in the data streams, and at the same time 
tests that a candidate GW signal appears in data from at 
least two interferometers [41| . Previously, we presented 
the results of a search for a GWB associated with the 
bright and nearby GRB 030329 [42]. Here, we present 
analysis methods which search for GWBs associated with 
GRBs not only on an individual- GRB basis to target loud 
GWBs, but also on a statistical basis. The statistical ap- 
proach is sensitive to the cumulative effect of any weak 
GW signals that may be present in the LIGO data. 

It is noted here that for the compact binary coalescence 
models of short-duration GRBs, a subset of the associ- 
ated inspiral waveforms are well modelled, and that a 
template-based search for inspiral GW signals associated 
with short-duration GRBs is currently being developed 
using LIGO data. 



II. LIGO S2, S3, AND S4 SCIENCE RUNS 

The LIGO interferometers (IFOs) have been described 
in detail elsewhere [43]. These detectors are kilometer- 
length Michelson interferometers with orthogonal Fabry- 
Perot arms, designed to detect impinging gravitational 
waves with frequencies ranging from ~ 40 Hz to sev- 
eral kilohcrtz. The interferometers' maximum sensitivity 
occurs near 100 Hz to 200 Hz. There are two LIGO ob- 
servatories: one located at Hanford, WA (LHO) and the 
other at Livingston, LA (LLO). There are two IFOs at 
LHO: one IFO with 4-km arms (HI) and the other with 
2-km arms (H2). The LLO observatory has one 4-km 
IFO (LI). The observatories are separated by a distance 
of 3000 km, corresponding to a time-of-ffight separation 
of 10 ms. 

Each IFO consists of mirrors at the ends of each arm 
which serve as test masses. Data from each IFO is in 
the form of a time series, digitized at 16384 samples/s, 
which records the differential length of the arms and 
which, when calibrated, measures the strain induced by 
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FIG. 1: Progression of LIGO sensitivities from S2 to S4 sci- 
ence runs. For each run, the corresponding curve gives the 
magnitude of the noise spectral density, in strain-equivalent 
units, for one of the IFOs during a representative time inter- 
val within the run. The solid curve gives the initial LIGO 
design sensitivity goal as given in LIGO's Science Require- 
ments Document (SRD). 

a gravitational-wave. The response of an IFO to a given 
strain is measured by injecting sinusoidal excitations with 
known amplitude into the test mass control systems and 
tracking the resulting signals at the measurement point 
throughout each run. The result is a measurement of the 
time-varying, frequency-dependent response function of 
each IFO. 

The LIGO S2 run was held from February to April 

2003 (59 days), the S3 run from October 2003 to January 

2004 (70 days), and the S4 run from February to March 

2005 (29 days). The sensitivity of the LIGO detectors 
improved significantly between the S2 and S4 runs, and 
approached the initial LIGO design sensitivity during the 
LIGO S4 run. The progression of the best LIGO sensi- 
tivity from the S2 to S4 runs is shown in Fig. Q] For each 
run, the corresponding curve in this plot gives the mag- 
nitude of the noise spectral density, in strain-equivalent 
units, for one of the IFOs during a representative time 
interval within the run. The solid curve gives the initial 
LIGO design sensitivity goal as given in LIGO's Science 
Requirements Document. Further, the duty factor of the 
three IFOs increased significantly from the S2 to S4 run. 
During the S2 run, the duty factors were 74%, 58%, and 
37% for the HI, H2, and LI IFOs, respectively, while 
during the S4 run, the duty factors were 80.5%, 81.4%, 
and 74.5%, respectively. 



III. THE GRB SAMPLE 

Compared to the 1990s, when BATSE was detecting 
GRBs, the period from 2001 to 2004, when LIGO had its 
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first three science runs, was a time of relatively low GRB 
detection rate. LIGO's S4 run coincided with a time 
when Swift had just started operating and was making 
its first GRB detections. There were 29 GRB triggers 
during the S2 run, 11 GRB triggers during S3, and 6 GRB 
triggers during S4. These GRB triggers were provided by 
the Third Inter-Planetary Network (IPN) [441 ]. Konus- 
Wind HI, HETE-2 , INTEGRAL, and Swift, and were 
distributed via the GRB Coordinates Network (GCN). 1 

Only LIGO data which were of science mode quality 
were analyzed. These science mode segments are data 
collected when the interferometers were in a stable, res- 
onant configuration. Additionally, data segments which 
were flagged as being of poor quality were not included 
in the analysis. For example, data segments which were 
known to have a high rate of seismic transients were ex- 
cluded from the analysis. After all the data quality cuts 
were made, there were 28 GRBs left to be analyzed for 
the S2 run, 7 GRBs for S3, and 4 GRBs for S4, for a to- 
tal of 39 GRB triggers. Of these, 22 GRBs had positions 
well-localized to within a few arcminutes, while 17 GRBs 
did not. These 17 GRBs were detected by either HETE 
or IPN. In the case of HETE, no position measurements 
were available while, in the case of IPN, the GRBs were 
not well-localized. Of the 39 GRBs, six had redshift mea- 
surements, four of which were at z > 1, and two fell in 
the short-duration category of bursts, i.e. had durations 
< 2 seconds. For this analysis, due to the small size of 
the sample, we did not attempt to differentiate the GRBs 
according to their observed properties. The use of a clas- 
sification scheme in a search can be done in the future 
with a larger GRB sample. 

Information about most of the GRBs was collected 
from the corresponding GCN circulars. The parameters 
that are relevant for this analysis are the GRB date and 
trigger time, and the right ascension and declination. For 
those HETE GRBs which did not have positions, infor- 
mation about the GRB trigger time was obtained from 
the HETE website. 2 A list of the GRBs analyzed and 
relevant information are given in Table |TJ 



IV. DATA ANALYSIS 

A. On-source and off-source data segments 

Since GRBs have well-measured detection times, the 
search for short-duration GW signals can be limited to 
time segments — called on-source segments here — sur- 
rounding the GRB trigger times. Limiting the search 
to encompass only these time segments significantly re- 
duces the number of search trials, compared to a search 
which makes use of data from the entire run. In case 



of a detection, such a reduction in trials translates to 
a larger significance for the detection compared to that 
which would result from an untriggered search. 

Making use of on-source segments also means that 
background estimation can proceed by using data 
stretches — called off-source segments here — which are 
outside the on-source segments, but which are still close 
enough in time to the on-source segments so that the off- 
source data are similar in character to, and representative 
of, the on-source data. 

In this analysis, the length of each on-source segment 
was chosen to be 180 seconds, with the first 120 seconds 
of the LIGO on-source data occurring before the GRB 
trigger time, and the last 60 seconds occurring after the 
trigger time. This window length is longer than the ex- 
pected time delay between a gravitational-wave signal 
and the onset of a GRB sig nal, which is of the order of 
several seconds [i^, |47], l4c| . but which in certain mod- 
els can be as large as 100 seconds [49(. The large search 
window also takes into account the uncertainty in the def- 
inition of the measured GRB trigger time, i.e. it takes 
into account the possibility that the trigger time used in 
the analysis occurred before or after the actual start of a 
gamma-ray burst signal. Many gamma-ray light curves 
show sub-threshold, precursor bursts which occur before 
the measured GRB trigger time, hence our choice of an 
asymmetric search window around the trigger time. 

For each GRB, a search for a GW signal was carried 
out using data from each pair of IFOs that was operating 
properly at that time. Additionally, LHO-LLO on-source 
pairs were analyzed only when GRBs had well-defined 
positions, since position information is necessary to cal- 
culate the LHO-LLO time-of-flight delay. After all the 
data quality cuts were made, there were 59 IFO-IFO on- 
source pairs that were analyzed. This number is larger 
than the number of GRB triggers because, for each GRB 
trigger, it was possible to have up to three IFO pairs pass 
the data quality cuts. There were 35 H1-H2 on-source 
pairs analyzed, 12 for Hl-Ll, and 12 for H2-L1. 

The software used in this analysis is available in the 
LIGO Scientific Collaboration's CVS archives with the 
tag multigrb_rl in MATAPPS. 3 



B. Data conditioning 

Before the crosscorrelation between two LIGO data 
streams was calculated, the time series data from each in- 
terferometer was conditioned. This consisted of whiten- 
ing, phase-correction, and bandpassing from 40 Hz to 
2000 Hz. The sampling rate was retained at 16384 sam- 
ples/s. Whitening was done to make sure the resulting 



1 http://gcn.gsfc.nasa.gov 

2 http://space.mit.edu/HETE 



http : / / www. lsc- group . phys . uwm . edu/cgi-bin/evs / viewc vs . cgi 
/ matapps /sre/ searches /burst/ multigrb 
/?cvsroot=lscsoft&sortby=rev#:dirlist 
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TABLE I: The GRB sample analyzed 
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"For GRBs with the same date, letters are appended to the date 
to distinguish the GRBs. 
6 UTC time of GRB trigger. 

C GPS time of GRB trigger (seconds since Oh 6 Jan 1980 UTC.) 
^Duration of gamma-ray burst. 
e Right Ascension of GRB. 
^Declination of GRB. 

"Polarization-averaged antenna factor for specified IFO site 
(cf. Eq. EJ. 

ft Time-of-fiight of GW signal between LHO and LLO. A positive 
value means that the signal arrived first at LLO; a negative value 
means that the signal arrived first at LHO. 

'Interferometers which were analyzed. 

H = 1.986. 

k z = 3.372. 
l z = 0.168. 

m z = 0.5915. 

n z = 1.44. 

°z = 3.24. 
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spectrum of the data was flat instead of being domi- 
nated by low- frequency or high-frequency components. 
The procedure consisted of using one-second data units 
to whiten the adjacent one-second data and, as a conse- 
quence, removed any non-stationarity in the data having 
a time scale larger than one second. The whitening pro- 
cedure also removed known lines. 

The response functions of the three LIGO interferome- 
ters to a given GW strain signal are not exactly the same. 
A GW signal impinging on the three interferometers will 
thus appear as having slightly different phases in the cor- 
responding time series data (even after correcting for the 
LHO-LLO time-of- flight delay) . Phase correction of the 
time series data was therefore done to remove the dif- 
ferences that can be attributed to the different response 
functions of the interferometers. The phase correction 
process made use of the measured, time-dependent, re- 
sponse functions of the interferometers. 



C. Measuring the crosscorrelation statistic 

The search method consisted of a simple "binned" 
search in which the 180-second conditioned on-source 
time-series for each IFO was divided into time inter- 
vals (or bins) and the crosscorrelation for each IFO-IFO 
time bin pair calculated. Crosscorrelation bins of lengths 
25 ms and 100 ms were used to target short-duration 
GW signals with durations of ~ 1 ms to ~ 100 ms. 
These crosscorrelation lengths were found, through sim- 
ulations, to provide sufficient coverage of the targeted 
short-duration GW signals. Using bins much shorter 
than 25 ms would considerably increase the trials in the 
search, and therefore decrease the significance of a can- 
didate GW event, while using bins much longer than 
100 ms would considerably diminish the crosscorrelation 
strength of signals in the two data streams due to the 
increased duration of noise. The crosscorrelation, cc, is 
defined as: 



5^[si(i) - fii][s 2 (i) -m] 



i=i 



\ 



£[«iCj)-mi] 5 

3=1 



(1) 



. £Mfc) - M2] 
\ fc=i 



where si and S2 are the two time series to be correlated, 
111 and H2 are the corresponding means, and m is the 
number of samples in the crosscorrelation, i.e. the cross- 
correlation integration length multiplied by the sampling 
rate of 16384 samples/s. The possible values of the nor- 
malized crosscorrelation range from —1 to +1. 

The bins were overlapped by half a bin width to 
avoid inefficiency in detecting signals occurring near a 
bin boundary. The crosscorrelation value was calculated 
for each IFO-IFO bin pair and, for each crosscorrela- 
tion bin length used, the largest crosscorrelation value 
in the case of an H1-H2 search — obtained within 



the 180-second search window was considered the most 
significant measurement for that search, for that cross- 
correlation bin length, for that IFO pair. In the case of 
an Hl-Ll or H2-L1 search, it was the largest absolute 
value of the crosscorrelations that was taken as the most 
significant measurement. This was done to take into ac- 
count the possibility that signals at LHO and LLO could 
be anti-correlated depending on the gravitational wave's 
(unknown) polarization. In the sections that follow, a 
reference to the "largest crosscorrelation" , in the case of 
an LHO-LLO analysis, will always mean the largest ab- 
solute value of crosscorrelations. 

For those GRBs which had well-defined positions, the 
position of the GRB in the sky at the time of the burst 
was used to calculate the GW signal's time-of-flight delay 
between the LHO and LLO observatories. Each LHO- 
LLO pair of 180-second on-source segments were shifted 
in time relative to each other by the corresponding time- 
of-flight amount before the crosscorrelations were calcu- 
lated. For those GRBs which were not well-localized, 
only H1-H2 on-source pairs were analyzed. For these 
GRBs, the maximum uncertainty in the LHO-LLO time 
delay is ±10 ms, which is of the same scale as the signal 
durations targeted by the analysis, and such a time offset 
between signals at the two interferometers would have a 
considerable effect on the measured crosscorrelation. 



D. Post-trials distributions 

To estimate the significance of the loudest event, i.e. 
the largest crosscorrelation, that was found in an on- 
source segment corresponding to a GRB and an IFO 
pair, we used off-source data within a few hours of the 
on-source data to measure the crosscorrelation distribu- 
tion of the noise. This distribution was obtained for each 
GRB, for each IFO pair, for each crosscorrelation length 
by applying the search (described in Sections IIVBI to 
IIV C[) on the off-source segments. The total length of the 
off-source region was about three hours surrounding the 
on-source segment. Each distribution was constructed by 
collecting the largest crosscorrelation (or largest absolute 
value of crosscorrelations, in the case of Hl-Ll and H2- 
Ll) from each 180-second segment of the off-source re- 
gion. This post-trials distribution takes into account the 
number of effective trials that was used in searching the 
on-source segment. 

To obtain enough statistics for each distribution, time 
shifts were performed such that the time series of each 
IFO was shifted by multiples of 180 seconds relative to 
the other IFO and two 180-second stretches from the two 
IFOs were paired at each shift, making sure that two 180- 
second time stretches were paired only once for each dis- 
tribution. The time shift procedure effectively increased 
the length of the off-source data to about 50 hours or 
more, typically. 

As an example, the post-trials distribution for GRB 
050318, for the H1-H2 IFO pair and for the 25-ms cross- 
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(cf. Section III), the applied data conditioning (cf. Sec- 
tion IVB), and the use of off-source data immediately sur- 
rounding the on-source data to estimate the background 
noise (cf. this section), which made it more likely that 
the background would properly reflect the rate of any 
correlated noise in the on-source data. 

The cumulative distribution of local probabilites re- 
sulting from the search of 59 on-source segment pairs is 
shown in Fig. [3] for the 25-ms crosscorrelation length, and 
in Fig. H] for the 100-ms crosscorrelation length. Also 
shown (bold dashed lines) is the expected distribution 
under a null hypothesis. There were no loud events that 
were not consistent with the expected distribution, and 
we therefore conclude that there was no loud GW signal 
associated with any single GRB in the sample. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

largest crosscorrelation 



STATISTICAL TESTS 



FIG. 2: (a) Example of a crosscorrelation post-trials distri- 
bution for the 25-ms crosscorrelation window, for the H1-H2 
IFO pair. Data from off-source segments for GRB 050318 was 
used, (b) Cumulative distribution of (a), normalized to the 
total number of entries in the distribution. Both distribu- 
tions with and without time shifts are shown, including the 
statistical errors. The arrow points to the largest crosscorre- 
lation found in the on-source segment for GRB 050318. In 
this example, the largest crosscorrelation of 0.36 has a local 
probability of 0.57. 



correlation length, is shown in Fig. O For comparison, 
the cumulative plot shows both the distribution obtained 
with time shifts, and the distribution obtained without 
employing time shifts. 

Each resulting post-trials distribution was used to es- 
timate the cumulative probability that the largest cross- 
correlation found in the corresponding on-source segment 
could be due to noise. This was done by determining 
what fraction of the distribution were at least as large as 
the loudest crosscorrelation found in the on-source seg- 
ment. For example, the significance of the loudest 25-ms 
crosscorrelation found in the H1-H2 on-source segment of 
GRB 050318, indicated by an arrow in Fig. [2Kb), can be 
estimated by using the plotted post-trials distribution. 
This probability will be referred to interchangeably in 
this paper as the post-trials, or local, probability of the 
on-source crosscorrelation statistic. This is also known 
in the literature as the false alarm probability. 

Since HI and H2 are colocated, environmental distur- 
bances can give rise to correlated transient noise in the 
two interferometers. The effect of these correlated envi- 
ronmental noise on an H1-H2 crosscorrelation were, how- 
ever, suppressed by: the judicious use of data quality cuts 



As mentioned earlier, GW signals from individual 
GRBs are likely to be weak in most cases due to the cos- 
mological distances involved. Therefore, besides search- 
ing for GW signals from each GRB, we also consider the 
detection of a GW signature associated with a sample 
of several GRBs. Such approaches, first proposed in the 
context of GWs in [5(3], have already been used fUHH to 
analyze resonant mass detector data using triggers from 
the BATSE and BeppoSAX missions. 

We use two different statistical methods to look for 
a GW signature associated with a sample of multiple 
GRBs. As one may expect, the statistical performance 
of a method will depend on the nature of the underlying 
source population distribution. The two different meth- 
ods presented here have complementary properties in this 
respect. The first statistical test presented, the binomial 
test, is most effective when several events contribute to 
the tail, i.e. the significant end, of the probability distri- 
bution of a sample. Moreover, it is also effective when 
there is a single significant event in the sample. The sec- 
ond test, the rank-sum test, is more effective at detecting 
the cumulative effect of weaker signals, but it is not very 
effective at detecting a few large events which fall on the 
tail of a probability distribution. 

Since the signal strengths targeted by these two meth- 
ods are slightly different, the resulting significances from 
the two methods can be different when there are real sig- 
nals present in the sample. If a detection is claimed and 
the more significant measurement from the two statistical 
tests is chosen, then the proper statistical treatment, in 
order to arrive at a final significance, would be to impose 
a penalty factor for using two statistical tests to search 
for the cumulative signal. 
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FIG. 3: Cumulative local probability distribution resulting 
from the search of 59 IFO-IFO on-source pairs using a 25-ms 
crosscorrelation length. The most significant excess is indi- 
cated by the arrow. The expected distribution under the null 
hypothesis is indicated by the bold, dashed line. The excess 
needed for a ~ 1% confidence in the null hypothesis is indi- 
cated by the solid line. The maximum excess indicated by this 
line is 15 events because only the 15 most significant events 
in the actual distribution are tested. 



Testing a probability distribution: The binomial 
test 



Under a null hypothesis, the distribution of local prob- 
abilities is expected to be uniformly distributed from to 
1. The measured distribution of local probabilities was 
tested to search for an excess which may have been due to 
the cumulative effect of weak GW signals. In particular, 
we searched the tail of the distribution, or the smallest 
probabilities found in the on-source searches, by using 
the binomial test. To test the tail of a probability distri- 
bution, one first makes a choice as to how many events, 
n, in the tail would be tested out of the total number 
of events, N, in the sample. In this analysis, there were 
59 IFO-IFO on-source pairs, and the upper 25% of the 
resulting probability sample, or the 15 most significant 
events, was tested. The probabilities of these n events 
are then sorted according to increasing value, i.e. de- 
creasing significance: pi,P2,P3, ■■-,Pi, ■■■,Pn- For each of 
these probabilities, pi, one calculates the cumulative bi- 
nomial probability, which is the probability for getting i 



FIG. 4: Similar to Fig. [3] but using a 100-ms crosscorrelation 
length. 



or more events at least as significant as Pi\ 

P>AP 1 ) = PiiPi) + Pi+i(Pi) + Pi+2(Pi) + ••• + Pfffa) 

(2) 

= 1 - [Po(Pi) + Pi(Pi) + P 2 (Pi) + .- + Pi-lfa)] 

(3) 

and where Pi(p) is the binomial probability for getting i 
successes in N trials: 



Piip) 



Nl 



il(N-i) 



-p\l-p) 



N- 



(4) 



Here, TV is the number of on-source searches, which is 59, 
and "success" means getting i events at least as signif- 
icant as p. Note that if there is one loud event in the 
sample, with p -C 1, then it follows from Eqs. [3] and [3] 
that the cumulative binomial probability is, 



P> 1 (p) = l-(l-p) N 
fa Np 



(5) 
(6) 



Thus, the binomial test is able to automatically handle 
the case of a single loud event in the distribution. 

After the cumulative binomial probability, P>i(j>i), has 
been calculated for each post-trials probability, p i: the 
smallest binomial probability in the set is identified. This 
smallest binomial probability will point to the most sig- 
nificant excess that was found in searching the tail of the 
probability distribution. 

The most significant excess that was found by the bi- 
nomial test in the tail of the distribution is indicated 
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by an arrow in Figs. [3] and [4j For the 25-ms dis- 
tribution, the smallest binomial probability found was 
P>9(P9 — 0.104) = 0.153. This means that the binomial 
test found that the most significant excess in the tail of 
the distribution consisted of nine events with local prob- 
abilities p < 0.104, and that the binomial probability for 
having nine or more events at least as significant as 0.104, 
given 59 trials, is 0.153. 

In the case of the 100-ms distribution, the smallest 
binomial probability found was P>g(pg — 0.112) = 0.207. 
This means that the binomial test found that the most 
significant excess in the tail of the distribution consisted 
of nine events with local probabilities p < 0.112, and that 
the binomial probability for having nine or more events 
at least as significant as 0.112, given 59 trials, is 0.207. 

Searching the tail of a post-trials probability distribu- 
tion for the most significant excess introduces additional 
trials to the search. We thus need to test the most sig- 
nificant excess found in the tail of each local probabil- 
ity distribution against the null hypothesis to properly 
establish its level of significance. The expected distribu- 
tion of the binomial probability statistic under the null 
hypothesis was obtained through simulations. The simu- 
lations consisted of randomly generating 59 numbers uni- 
formly distributed from to 1 to simulate 59 post-trials 
probabilities under the null hypothesis. Then the same 
binomial test that was applied to the actual post-trials 
probability distribution was applied to this distribution 
of random events to search for the most significant ex- 
cess in the 15 most significant events in the tail. This 
was repeated a million times, and the binomial probabil- 
ity of the most significant excess found in each trial was 
collected. The resulting distribution of binomial prob- 
abilities under the null hypothesis, in effect, takes into 
account the number of trials used in searching the tail of 
the post-trials distribution. 

Results of these simulations show that, under the null 
hypothesis, the probability for getting a measurement at 
least as significant as 0.153 that was found in the 25-ms 
search is 0.48. In other words, under the null hypothesis, 
1 in 2.1 sets of 59 on-source searches will result in a most 
significant excess with a binomial probability at least as 
significant as 0.153. This quantifies the conclusion that 
the result of the 25-ms search is consistent with the null 
hypothesis. 

Similarly, we find that, under the null hypothesis, the 
probability for getting a measurement at least as signifi- 
cant as 0.207 that was found in the 100-ms search is 0.58. 
In other words, under the null hypothesis, 1 in 1.7 sets 
of 59 on-source searches will result in a most significant 
excess with a binomial probability at least as significant 
as 0.207. And, as with the 25-ms result, this level of sig- 
nificance for the 100-ms search result is consistent with 
the null hypothesis. 

Also shown in Figs. [3] and [4] is a curve indicating the 
excess needed for a ~ 1% confidence in the null hypoth- 
esis. At each local probability, the curve gives the cumu- 
lative number of events needed to obtain a ~ 1% final 



probability under the null hypothesis, given 59 on-source 
pairs. 



B. Maximum likelihood ratio based tests 

A maximum likelihood ratio test [53| for detecting a 
GW signature associated with a sample of multiple trig- 
gers was derived in (54|. (It was shown there that [5(| 
is a special case of the maximum likelihood ratio ap- 
proach.) The method proposed in [54[ cannot be ap- 
plied directly to the entire GRB sample described above 
since the largest crosscorrelation values were obtained in 
diffe rent w ays for H1-H2 and Hl-Ll (H2-L1) (cf. Sec- 
tion [IVC|. In the following, we will only use the largest 
crosscorrelations from H1-H2 on-source segments. This 
reduces the total number of GRB on-source segments 
used in this test to 35. 

Let the largest crosscorrelation from the i th GRB on- 
source segment be denoted as cc max ,i- If we do not use 
any prior probability distribution for the properties of 
GW signals associated with GRBs, the maximum likeli- 
hood ratio detection statistic is simply the average of the 
largest crosscorrelation values from the GRB set, 

where Nqrb is the number of H1-H2 GRB on-source seg- 
ments used. We call \ the sum-max statistic. 

To build in robustness against instrumental noise arte- 
facts, such as short duration transients, we replace the 
sum-max statistic, which was derived for the ideal case 
of Gaussian and stationary noise, by a non-parametric 
counterpart. The on-source and off-source largest cross- 
correlation values are pooled into two separate sets and 
the Wilcoxon rank-sum test [55| is used for the null hy- 
pothesis that the two sets of samples were drawn from 
the same underlying true distribution. 

The cumulative distribution of the on-source and off- 
source largest crosscorrelations from the 100-ms search 
are shown in Fig. O Application of the rank-sum test 
shows that the significance of the null hypothesis is 0.64. 
This implies that one out of 1.6 trials can show a false 
positive detection at this significance threshold. Assum- 
ing that GRB triggers occur at a rate of one per day, 
one year of observation would contain approximately 10 
collections of 35 GRBs. In order to achieve a low false 
detection probability, we would require a much lower sig- 
nificance, such as < 0.01, in order to reject the null hy- 
pothesis. 

As a further check, we also compute the empirical sig- 
nificance of the on-source value of x with respect to the 
set of off-source x values. Values of the off-source % 
statistic were calculated by pooling the largest crosscorre- 
lations from the H1-H2 off-source segments, then dividing 
this pool into subsets, each of which had -/Vqrb number 
of elements. For each of these subsets, the x statistic 
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FIG. 5: The cumulative distributions of the on-source (solid 
black with + marker) and off-source (solid gray) largest Hl- 
H2 crosscorrelations from the 100-ms search. The vertical 
lines denote the locations of the medians of the off-source 
(gray) and on-source (black, dashed) samples. 
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FIG. 6: Plot of the cumulative distribution function of the 
off-source values of the sum-max statistic \- The dashed line 
shows the best fit normal distribution. The solid horizontal 
and vertical lines indicate the location of the on-source values 
of Xon and the corresponding cumulative probability. The 
on-source value was Xon = 0.1753, which yields a cumulative 
probability of 0.64 or an empirical significance of 1 - 0.64 = 
0.36. 



was calculated using Eq. [7J The empirical significance is 
defined as the fraction of off-source x values greater than 
or equal to the on-source x value. The empirical sig- 
nificance has more scatter than a significance computed 
with a known null distribution since we only have a finite 
number of off-source values. However, the number of off- 



source values in this analysis is large enough that we can 
ignore the associated error. 

The mean and standard deviation of the off-source 
sum-max sample were 0.1744 and 0.0029, repectively. 
The on-source value of sum-max was % on = 0.1753. Fig- 
ure [H] shows the distribution of the off-source values of 
the test statistic. The empirical significance of the null 
hypothesis is 0.36. Following the discussion above, it is 
clear that this result is consistent with no detection. 



VI. SINGLE-GRB LIMITS 

Simulations were done to estimate the sensitivity of the 
search method to incident GW burst signals. This pro- 
cess was limited by the fact that the theoretical wave- 
forms of the GW burst signals associated with GRBs 
were not known. Other unknown quantities were: the 
polarization of the waves, the orientation of the source 
relative to the observer, and the redshifts of most of the 
GRBs. Conscious of these limitations, we proceed to set 
upper limits on the root-sum-square amplitude (h iaa ) of 
GW burst signals incident on the interferometers during 
the on-source times by using simulated waveforms with 
burst-like characteristics, adding these waveforms to the 
raw IFO data streams, and measuring the resulting cross- 
correlations. 

The antenna response of an IFO to incident, inde- 
pendent gravitational- wave strains, h+(t) and h x (t), de- 
pends on the relative position of the source in the sky 
and the polarization of the wave [56| : 

h(t) = F+(6, 0, VOM*) + (0, 0, vofc* (t) , (8) 

where (6, 4>) is the position of the source relative to the 
IFO's zenith and x-arm, respectively; ip is the polar- 
ization angle of the gravitational-wave; and F + (9, (j>,i/}) 
F x (8, (j>, are the corresponding "plus" and "cross" an- 
tenna factors. For most of the GRBs analyzed, the posi- 
tion, (#,</>), was known. The polarization angle, tp, how- 
ever, was an unknown parameter for all of the GRBs. 
Since the antenna factor is used in the simulations, upper 
limits were not set for GRBs which did not have well- 
defined positions. The polarization-averaged antenna 
factor is defined as: 



n 



F 2 



<*?>W(*s>*- ( fl ) 



We used sine-gaussians as the simulated waveforms for 
h + (t) and cosine-gaussians for h x (t) in Eq. [5) 



h + (t) = h+ o sin(27r/o£) exp 
h x (t) = /ix,ocos(27r/ i)exp 



2Q 2 

-(27T/ t) 2 

2Q 2 



(10) 
(11) 



where /o is the central frequency of the sine-gaussian and 
cosine-gaussian, h + and /i x ,o are the amplitude parame- 
ters of the + and x polarization signals, respectively, and 
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Q is a dimensionless constant which represents roughly 
the number of cycles with which the waveform oscillates 
with more than half of the peak amplitude. The root- 
sum-squared (rss) amplitude of h + (t) and h x (t) is related 
to these parameters via: 



\h+{t)\> dt « h+,o 



\h x (t)\ 2 dt « h xfi 



for Q > 3 , (12) 



for Q > 3 . (13) 



Using these waveforms for h + (t) and h x (t), we simu- 
lated circularly polarized GW waves by setting the sine- 
gaussian and cosine-gaussian amplitudes equal to each 
other, /i+.o = h x ,o = ^o- To simulate linearly polarized 
waves, we set /i x ,o = 0. In the discussion that follows, 
we define the h TSS of a simulated waveform as: 



(M*)l a + IM*)P)# 



(14) 



Since the polarization angle, tp, was not known for any 
GRB, a random polarization angle from to 360 degrees 
was generated for each simulated waveform event. In 
the case of LHO-LLO simulations, the source position- 
dependent difference in the polarization angles at LHO 
and LLO — due to the non-aligned detector arms — was 
taken into account. Finally, before adding the attenuated 
waveform given by Eq.[8]into an IFO's raw data stream, it 
was first calibrated using the measured response function 
ofthelFO. 

Following the procedure outlined above, simulated 
sine-gaussians with different frequencies and h rss values 
were added to each IFO on-source data at known times. 
Randomness in the injection times of the order of the 
crosscorrelation length was introduced to simulate the 
fact that the relative time of the GW signal within the 
180-second search window was not known. Data with in- 
jected signals were then conditioned using the procedure 
outlined in Section IIVBI The search was then applied 
to the data near the injection times — not to the entire 
180-second on-source segment — to find the largest cross- 
correlations around the injection times. This simulation 
procedure resulted in the determination of the probabil- 
ity density, p(cc\h TSB ), for measuring a crosscorrelation, 
cc, corresponding to a signal injected in an on-source seg- 
ment with a certain h rss value. 

The method used to set upper limits on h rss follows 
the standard recipe for setting frequentist upper limits 
[5^ . If p(cc\h rss ) is the probability density for measuring 
a crosscorrelation, cc, in an on-source segment given a 
signal with a certain h TSS value, then the 90% upper limit 
curve can be constructed from the set (ccgo, h ISS ), such 
that, 



0.90 



p(cc\h laa ) d{cc) 



(15) 
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FIG. 7: Examples of upper limit curves that were used to 
set upper limits on h ISS using linear and circular polariza- 
tions. These were the upper limit curves for the H1-H2 IFO 
pair, for GRB 050306, using sine-gaussians with Q = 8.9 and 
fo — 150 Hz. The shaded regions indicate the total la uncer- 
tainty in the h raa value. 



Examples of upper limit curves obtained through this 
procedure are shown in Fig. with one curve correspond- 
ing to linear polarization, and the other curve corre- 
sponding to circular polarization. These curves were ob- 
tained using the H1-H2 on-source data for GRB 050306; 
150-Hz, Q = 8.9 sine-gaussians; and a 25-ms crosscorrela- 
tion length. Each curve shows the h rss value of the simu- 
lated waveform versus ccgo , the crosscorrelation value at 
which 90% of the measured crosscorrelation values were 
larger (see Eq. [T5")) . The data was fitted with a four- 
parameter sigmoid function, 



cc 90 =Pi + 



I -Pi 



p 4 (l + exp[-p 2 (log 10 (/i rss ) -P3)}) 



, (16) 



where parameter pi defined the asymptote of CC90 at 
small values of h ISS , p^ tracked the asymptote of ccgo 
at large values of h rss (i.e. P4 « 1/asymptote), p% was 
the value of h ISS which gave a mid-range value of ccgo, 
and P2 defined the slope of the curve. The largest cross- 
correlation found in the on-source segment is also shown 
in Fig. [7] (vertical dashed line). The 90% h xss upper limit, 
before uncertainties, was found by evaluating the upper 
limit curve, which is the inverse of Eq. [TBJ at the largest 
on-source crosscorrelation value found in the search. 

The curves in Fig. [7] also show the estimated total la 
uncertainty in the measurement of the h rss values. The 
uncertainty in the h ISS values comes from measured ran- 
dom and sytematic errors in the calibration parameters 
that were used to calibrate the simulated waveforms, and 
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FIG. 8: Progression of h ISS upper limits from the S2 to S4 
LIGO runs for linearly polarized sine-gaussian waveforms; 
25-ms crosscorrelation. 
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TABLE II: Best 90% h las upper limits resulting from a search 
of GW signals from GRBs occurring during the three LIGO 
runs; 25-ms crosscorrelation analysis (Hz _1//2 ) 



Run 



ftrss,90 

(circular) 



/lraa,90 

(linear) 



S4 
S3 
S2 



1.1 x 10 
8.5 x 10 

8.2 x 10 



" 21 (050306) 
■ 21 (031109A) 
' 21 (030414) 



3.6 x 10~ 21 (050223) 150 Hz 
2.9 x 10 - 20 (031109A) 250 Hz 
3.1 x 10~ 20 (030329B) 250 Hz 



It can be seen that the upper limits for the two cross- 
correlation lengths do not differ much for the waveforms 
that were used. The upper limits for /o = 250 Hz and 
25-ms crosscorrelation length are plotted in Figs. [8] and 
[9] for linear and circular polarizations, respectively. The 
improvement in sensitivity from the S2 to S4 runs can be 
seen in these plots. The best upper limits from the three 
science runs are given in Table ITT1 From the S2 to the S4 
run, there was an improvement in sensitivity by about 
an order of magnitude. 

It can also be seen from Figs.[8]and[9]that, for most of 
the GRB source positions, the circular polarization limits 
are better than the linear polarization limits by about a 
factor of 3.5. This is always true in the case of H1-H2 
upper limits since waveforms at the two co-aligned LHO 
IFOs were always in phase (after calibrations). For LHO- 
LLO upper limits, there were two cases, GRB 030217 and 
030323a, in which the positions of the GRBs relative to 
the IFOs were such that circularly polarized waveforms at 
LHO and LLO were sufficiently out of phase so that up- 
per limits for circular polarization were not determinable 
for those GRBs. 



FIG. 9: Progression of h ISS upper limits from the S2 to S4 
LIGO runs for circularly polarized sine-gaussian waveforms; 
25-ms crosscorrelation. 



VII. CONSTRAINING GRB POPULATION 
MODELS 



also from the statistical errors which come from the sim- 
ulation procedure. Depending on which science run and 
IFO pair is being considered, the total ler uncertainty 
from all these sources ranged from ~ 10% to ~ 13%. 
However, for GRB 030217 and GRB 030226, the total 
uncertainty was about ~ 22% for the H1-H2 and Hl- 
Ll IFO pairs, due to larger calibration errors during the 
times of those GRBs. The final 90% h ISS upper limits 
were obtained by adding the corresponding total 1.28a 
uncertainties to the values obtained from the upper limit 
curves. 

The upper limits resulting from the use of Q = 8.9 sine- 
gaussians and a 25-ms crosscorrelation length, for GRBs 
with well-localized positions, are listed in Tables IIIII to 
[V]for linearly polarized waveforms, and in Tables IVT1 to 
IVIIII for circularly polarized waveforms. Corresponding 
limits from the use of a 100-ms crosscorrelation length 
are listed in Tables [IX] to EB and in Tables ED to EV] 



The approach of combining multiple GRBs to look for 
a GW signature associated with a sample of GRBs was 
described in Section [V] Having established that the null 
hypothesis is quite significant, i.e., that we cannot claim 
the detection of an association between GWs and the 
GRB population at a high enough confidence, we turn to 
setting constraints on the parameters of GRB population 
models. The method is summarized below and described 
in detail in [58j |. 

For a pair of detectors, it can be shown that only three 
scalar parameters associated with a GW signal are suffi- 
cient to determine the distribution of largest crosscorre- 
lations. The parameters are the matched filtering signal- 
to-noise ratios (SNRs) of the strain signals in individual 
detectors and the angle between the two strain signal vec- 
tors (as defined by the Euclidean inner product). In the 
following, a source population model is the joint proba- 
bility distribution function of these three parameters. 

Our approach to putting constraints on source popula- 
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tion models follows the standard frequentist upper limit 
procedure (cf. Section IVij) . In this case, let P(xl^source) 
be the marginal cumulative probability density function 
of the sum-max statistic, given the population model 

^source, and let Xa be SUCh that P{x a | -^source) = 1 - a, 

where < a < 1, and 1 — a is the desired confidence level. 
If the observed value of \ is greater than or equal to % Q , 
the corresponding model Z sourcc is accepted. It is rejected 
when x < Xa- To obtain the marginal distribution of 
we first construct its conditional distribution for a set 
of -^Vgrb values for the scalar parameters above, where 
Agrb is the number of H1-H2 GRB on-source segments. 
The marginal distribution of x for a given source popula- 
tion model can then be estimated by randomly drawing 
values of the scalar parameters followed by drawing x 
from the corresponding conditional distribution. 

Since we use only the H1-H2 pair, which are perfectly 
aligned, the angle between the strain responses is zero. 
Further, for narrowband signals, the SNR values for HI 
and H2 can be related by the measurable ratio of their 
calibrated noise power spectral densities (PSDs). Hence, 
only one parameter, which we chose to be the SNR, p, of 
the signal in HI, is required. Thus, the source population 
model, Z source , is simply the univariate distribution of p. 
An additional point that needs to be accounted for is the 
variation in the sensitivities of HI and H2, both within 
the runs as well as the significant improvements from 
one run to the next. This is done by fixing a fiducial 
noise PSD, S {a) (f), and approximating the PSD of HI 
for each GRB as simply a scaled version of it. We set 
the fiducial noise PSD to the one corresponding to the 
initial LIGO design sensitivity for the 4-km IFOs 4 and 
compute the scale factor at a fixed frequency of 200 Hz, 
which was approximately where most PSDs had their 
minimum during the S2, S3, and S4 runs. 

We use the theoretical prediction of the observed red- 
shift distribution of GRBs given in [5^] to construct 
■^source (prediction for the scenario of star formation via 
atomic hydrogen cooling) . An alternative is to simply use 
the measured redshift distribution but [3l|, [6(| show that 
there is a significant selection bias that affects the mea- 
sured redshifts for Swift and non-Swift GRBs, both of 
which are used in our analysis. The model in [591 ] is valid 
for long-duration GRBs, which are expected to trace the 
massive star formation rate of the Universe. We fit a 
piecewise parabolic curve (with 3 pieces) to figure 1 of 
[59( and then use the same subsequent calculational steps 
given in to obtain the redshift distribution for a flux- 
limited detector such as Swift. Fitting the star formation 
rate with a smooth curve allows us to extend the red- 
shift distribution reliably to very small values of the red- 
shift. Having obtained the distribution, we directly draw 
random values of the redshift, z, from it. Each redshift 
value is then converted to the corresponding luminosity 



4 http://www.ligo.caltech.edu/~lazz/distribution /LSC_Data/ 



distance D (corresponding to a Friedmann-Robertson- 
Walker cosmological model with f2 m = 0.3, J7a = 0.7 
and H Q = 72 km sec -1 Mpc -1 ). 

A simple model is used for the GW emission from 
GRBs. We assume that GRBs are standard candles in 
GW that emit a fixed amount of energy, £qw, isotrop- 
ically with similar amounts of radiation in the two un- 
correlated polarizations + and x . Further, neglecting 
the effect of redshift on the signal spectrum, we assume 
that the spectra of the received signals h + and h x are 
centered at a fixed frequency of f a in a band that is suffi- 
ciently narrow such that the noise power spectral density 
is approximately constant over it. In this case, the SNR 
is given by 

p ^ 2F -7mm' (17) 

where we have expressed the SNR with respect to the 
fiducial noise PSD. Since the emission is isotropic, the 
energy emitted in gravitational waves is (cf. Section 
VIIIA), 

E GW ~ —q - 1 — ~ foKss ■ ( 18 ) 

To convert the luminosity distance, D, for a given GRB 
into SNR p, we use the normalization 

p = V2F avcP0 ^(i±^) 3/2 , (19) 
D \l + z J 

where Dq is chosen to be the most probable luminos- 
ity distance, at the corresponding redshift zq, and po 
is the observed SNR for a GRB that occured at Dq 
with an optimal sky location and the above properties 
for h + , h y and -Egw- The redshift distribution pre- 
dicted in [59|] for Swift has a peak at z — 1.8, which 
yields Dq = 13.286 Gpc. The acceptance-rejection rule 
above simply becomes an upper limit on the value of pq. 
Note that, because of the scaling of noise PSDs discussed 
above, po should be understood as the SNR of the strain 
response (for a GRB directly above the detector) that op- 
erates at design sensitivity. For GRBs that do not have 
direction information, random values for F ave are drawn 
from a uniform distribution on the celestial sphere. 

Finally, in terms of the upper limit, p U ppcr, obtained 
on poj we get an upper limit on Eqw, 

E GW < ^-^Lf SW(f )p* . (20) 



For z = 1.8, f = 200 Hz, and y/S^(f ) = 2.98 x 
10- 23 Hz~ 1/2 , we get E GW < 8.43 x 10 55 p 2 ppor ergs 
(= 47.3p 2 r Af c 2 ). 

Figure [10] shows the 90% upper limit confidence belt 
for pq. The on-source value of sum-max was x — 0.1753 
for the S2, S3, S4 GRB sample. Hence, po < 35.5 and 
E G w < 5.96 x 10 4 M Q c 2 . This limit is too high to be 
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FIG. 10: Upper limit confidence belts at 90% confidence level 
on po, the SNR at the most probable redshift for Swift GRBs 
given in [H . The solid line is the curve for all S2, S3, S4 GRBs 
that were used in the H1-H2 search (on-source x = 0.1753). 
The dashed line is the curve when only the last 10 GRBs 
from the above set are selected (on-source \ — 0.1702). The 
line with filled circles is for a hypothetical scenario with 35 
GRBs, all with an optimal sky location, and two identical 
and constant sensitivity detectors. The shifting of the curves 
horizontally is due to the change in the variance of \ as the 
number of GRBs is changed. For each value of po, 10,000 
values of x were drawn from its marginal distribution. 

of any astrophysical importance. However, as discussed 
later, future analyses may be able to improve by orders 
of magnitude on this result. 

Since the detectors during the S2 run were much less 
sensitive than S4, one may expect that dropping the S2 
GRBs from the analysis can improve the upper limit. 
Figure [10] shows the 90% level upper limit belt ob- 
tained for the case when only the last 10 GRBs, span- 
ning the whole of S4 and part of S3, were retained in 
the analysis. The corresponding value of x = 0.1702 
yields an upper limit of 24.6 on po- Thus, we obtain 
Egw < 2.86 x 10 4 Mqc 2 . This shows, as expected, that 
making judicious cuts on the sample of GRBs can lead 
to improvements in upper limits. The upper limit can 
probably be improved further by retaining only the S4 
GRBs, but for a small number of GRBs the distribution 
of x used is not valid and a more accurate calculation 
has to be done. In Fig. [TUl we also show the upper limit 
confidence belt for a hypothetical scenario that is likely 
for the ongoing S5 run: a sample size of about 35 GRBs 
with the HI and LI detectors operating at design sensi- 
tivity. The implications of this curve are discussed in the 
next section. 

The confidence belt construction outlined in this sec- 
tion is for illustrative purposes only. In particular, we 
have not taken into account factors such as (i) changing 
noise spectral shapes, (ii) red-shifting of the standard 



candle (K-correction) and possible systematic errors as- 
sociated with the population model used. A more com- 
prehensive investigation is planned for the S5 data. 

VIII. DISCUSSION 

This search is not very restrictive with respect to mod- 
els for astrophysical systems which give rise to GRBs. 
The main assumption we have made is that the GW emis- 
sion is limited in duration — we sum over periods of up 
to 100 ms, which is much greater than the characteris- 
tic times expected for GW burst emission in most GRB 
models. Given the LIGO sensitivity at the time of this 
search, it is not surprising that our experimental lim- 
its in this search do not place significant restrictions on 
the astrophysical models at present. However, given the 
rapid development of the field, it is not precluded that 
the limits presented here will provide guidance to GRB 
astrophysics in the near future. In any case, it is use- 
ful to get a sense for the interplay between the measured 
gravitational- wave strain limits for individual GRBs from 
Section [VI] and astrophysical models. So in this section 
we provide some astrophysical context to our experimen- 
tal limits. We emphasize that the estimates given below 
are for illustration, and are not to be construed as mea- 
sured astrophysical limits. 

The local gravitational-wave energy flux in the two in- 
dependent polarizations, h + (t) and h x (t), is l6lT] 

which can be integrated over the duration of a burst of 
gravitational radiation and over a closed surface to relate 
the strains evaluated on the surface to the total intrinsic 
energy associated with a source within this volume. For 
a source at the center of a sphere of radius r at negligible 
redshift, then dA = r 2 dfl, as usual. 

Since many of the GRBs in the sample are found 
to have significant redshifts, it is useful to generalize 
the above to cosmological distances. In this case, we 
can use the luminosity distance, D, which by definition 
relates the intrinsic luminosity, C, of an isotropically 
emitting source to the energy flux T at a detector by 
C = T(A-kD 2 ). For a non-isotropic emitter, we replace 
the 47T by an integration over solid angle. We note that 
T is by definition the left hand side of Eq. [2TI and the in- 
trinsic luminosity is C = dE e /dt e . To integrate this over 
the signal duration at the detector, we use dt = (l+z)dt e . 
Hence, the energy emitted in gravitational radiation is, 

E e = J dQ J T dt 

= w*i&J" a l{ h + + h ^ il (22 » 

For negligible redshifts, D = r. We note that D — D(z) 
is itself a function of the redshift, and in general depends 
on the cosmological model. 
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If the signal power at the detectors is dominated by a 
frequency f Q , as is the case for the sine-gaussian wave- 
forms introduced earlier, then Eq. [22] can be written in 
the approximate form 



We then find an analogous expression to Eq. 



E, 



IT C 



D 2 



4 G 1 + z 



dt 



(23) 



which allows a direct relation between E e and the observ- 
able h xss (see Eqs. IT2l and [T3|) . For sine-gaussian wave- 
forms, the approximation is quite good for Q > 3; the 
error is approximately 1/(1 + 2Q 2 ). We will assume here 
that the simulated waveforms are effectively local to the 
detectors. Specifically, the frequency f a is the measured 
frequency (which is related to the source frequency f e by 
fo = /e/(l + z)). Of course, some fraction of the source 
power might be shifted in or out of the sensitive LIGO 
band in frequency or expanded in time beyond our inte- 
gration time. We ignore any such effect here. 



A. Case I: Isotropic emission 

First, we consider a simple, but unphysical, example 
where the radiation is emitted isotropically, with equal 
power in the (uncorrelated) + and x polarizations. In 
this case, Eq. [23] becomes 



7r 2 r 3 D 2 
iso ~ G l + z /oftrss 



(24) 



Then for a LIGO sensitivity for some waveform repre- 
sented by h ISS , we might hope to be sensitive to a dis- 
tance 



D S3 70 Mpc 



100 Hz \ / E s 



U ) \M Q c 2 



1/2 



x (1 + z) 1 / 2 



(25) 



for an isotropic source which emits gravitational- wave en- 
ergy E iBO (in units of solar rest energy) at detected fre- 
quency f a . 



B. Case II: Long-duration GRBs 

For long-duration GRBs, we consider the scenario 
where such events are associated with a core collapse, 
perhaps involving a very massive progenitor [62l | . Gravi- 
tational wave emission has been simulated for supernova 
core collapse models for relatively light (~ 10M Q ) pro- 
genitors, for example, in Refs. [3g, [37(. These models 
invoke axisymmetry, with linearly polarized strain that 
is proportional to sin 2 i, where l is the angle with respect 
to the symmetry axis. 

Integrating over the full solid angle, Eq . |2"31 becomes 



E„ 
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2 
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(26) 



D S3 1 Mpc 



/ 100 Hz 
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As described earlier, our experimental limits correctly ac- 
count for the antenna pattern associated with each GRB. 
Hence, no additional factors are required in the equation 
above if one were to use values from the tables of results. 
However, if one wished, for example, to apply a theoreti- 
cal h TSS to a particular GRB, the antenna factors for each 
GRB are given in Table [J 

Core collapse simulations indicate that most of the 
gravitational radiation is emitted from the core bounce, 
and that E sn should be at most 10~ 7 Mc 2 (37[, or even 
smaller (3(| . For the very massive progenitors often asso- 
ciated with long-duration GRBs, the collapse process is 
uncertain. Whether there is a bounce at all, or simply a 
direct collapse to a black hole, depends [63| on the mass, 
metallicity, and angular momentum of the progenitor. In 
any case, there is no reason to believe that the efficiency 
for converting the collapse into gravitational radiation 
increases with the progenitor mass. 

In fact, the situation for GW detection in this sce- 
nario is especially unpromising. It is natural to align the 
symmetry axis of the (rotating) core collapse with the 
direction of the gamma-ray beam. Hence, i — would 
be along the line of sight to the detectors. For a typical 
gamma-ray beaming angle of half- width ~ 10°, then at 
best, where the detectors are at the edge of the beam, 
this would give a suppression factor of ~ 30. Finally, we 
note that long-duration GRBs are distant objects, with 
mean observed redshift of S3 2. 4. 5 Given their rcdshift 
distribution, the simulations to date indicate that detec- 
tion of long-duration GRBs is unlikely if core bounce is 
the dominant radiation mechanism. 

However, core collapse can potentially drive other 
mechanisms more favorable for gravitational radiation 
detection. In particular, bar mode instabilities are po- 
tentially very efficient radiators and do not suffer from 
the unfavorable alignment noted above for axisymmet- 
ric core bounces. Similarly, core fragmentation during 
collapse can lead to GW radiation from the inspiraling 
fragments. Reference |63| has examined these possibili- 
ties, and while the likelihood of bar instabilities or core 
fragmentation, along with their detailed properties, is un- 
certain, the resulting gravitational radiation is plausibly 
detectable for a nearby GRB. In such cases, Eqs. [2^1 and 
[3D]might be more appropriate descriptions of the radiated 
energy and distance to which we can detect the source. 

The nearest known GRB to date is long-duration burst 
GRB 980425 at D = 35 Mpc. From Eq. [23 LIGO de- 
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tection at 35 Mpc by the method described in this paper 
would require an efficiency of at least -E sn /M c 2 « 10% 
for a IMq system, much larger than the efficiency ex- 
pected from conventional core collapse, but perhaps not 
unreasonable in case of bar instabilities or core fragmen- 
tation. Unfortunately, the data considered here did not 
include any such nearby events. For example, during the 
(most sensitive) S4 run, the GRB sample consisted of 
only 4 events, all long-duration GRBs. The most nearby 
of these with a measured redshift was GRB 050223 
(z = 0.5915) at J) k 3.5 Gpc. Assuming linear polar- 
ization, we can obtain an estimate for sensitivity from 
the 90% upper limit for GRB 050223 from Table [TTJ This 
gives for E sn the value 1.6 x 10 4 M Q c 2 . This is in fact 
very close to the source luminosity maximum of c 5 /G 
which gives 2 x 10 4 M Q c 2 if sustained for 100 ms. 
The larger sample of GRBs in future runs will hopefully 
include some long-duration GRBs at smaller redshift. 



C. Case III: Short-duration GRBs 

Short-duration GRBs, to the extent that the popula- 
tion is associated with the merger of compact binary sys- 
tems, offer several potentially interesting characteristics. 
First, such mergers are found to be relatively efficient ra- 
diators of gravitational radiation. Second, the emission 
pattern is not expected to be problematic. Moreover, the 
measured redshifts to date indicate a significant num- 
ber of relatively low-z GRBs. (The average redshift was 
w 0.4 for the 2005 sample of 5 events.) The mergers may 
include formation of a hypermassive neutron star [65[ or 
a black hole with associated ringdown [6||. Finally, the 
merger GW emission, which is best suited to the method- 
ology described in this paper, would be preceded by an 
inspiral which is potentially detectable by a sensitive, in- 
dependent LIGO search based on matching inspiral wave- 
form templates [67] . However, we have verified that the 
present search, while not as sensitive to inspirals as the 
dedicated waveform template-based search, can readily 
detect inspiral emission when there is sufficient signal to 
background in individual 25-ms or 100-ms bins. In this 
case, the maximum crosscorrelation occurs when the fre- 
quency of the inspiral radiation passes through the 100- 
300 Hz range, where the detector sensitivity is best (see 
Fig. [J). Therefore, while this search is uniquely sensitive 
to the higher frequency, short-duration, poorly modeled 
gravitational waves from the merger phase, it also pro- 
vides independent information on the inspiral phase. Re- 
cent estimates [H], HH] place the chance for detection of 
a BH-NS merger at up to ~ 30% for a year of simultane- 
ous LIGO and Swift operation, and ~ 10% for a NS-NS 
merger. Here, we provide an estimate for a contrived, 
but physically motivated, model. 

We suppose that the gravitational-wave emission pat- 
tern for the merger follows that of the inspiral, that is 

h+ = h a f(t) -(1 + cos 2 i) ; h x = h a g(t) cost (28) 



where i is measured with respect to the axis orthogo- 
nal to the plane of the inspiral orbit. The functions 
f(t) and g(t) are orthogonal functions, for example f(t) 
could be the sine-gaussian form discussed earlier, while 
g(t) is a cosine-gaussian; h represents a constant am- 
plitude. While the degree of gamma-ray beaming for 
short-duration GRBs is still uncertain, we suppose that 
the gamma rays are preferentially emitted along the an- 
gular momentum axis of the merger system. If the Earth 
is near the center of the gamma-ray beam, then i = 
is along the line of sight between detector and source, 
which is a maximum of the assumed emission pattern, 
and the radiation will be circularly polarized. Returning 
to Eq.l23"land integrating over the full solid angle, we find 
in this case (with u = 0): 



27r 2 3 r>2 
rr ~ £__ __ f2 h 2 

Emerge ~ 5 G 1 + z Jo n rs 

Rewriting this for D, as before, gives 



(29) 



D « 44 Mpc 



/ 250Hz ^ (E mcrgc ^ 1/2 



10" 21 Hz 



-1/2 



(1 + *) 



1/2 



(30) 



The comments below Eq. [27] concerning antenna factors 
also apply here. 

There has been substantial recent progress in calcu- 
lations of gravitational radiation production in various 
types of mergers. Numerical simulations of NS-NS merg- 
ers give [g6|, [H, H^] typical values of the radiated energy 
of about 0.5-1% of the total mass, or E m 0.01M Q c 2 . 
These simulations indicate that the frequency spectrum 
can be broad, ranging from a few hundred Hz to ~ 2 kHz. 
Perhaps the most interestin g ca se is BH-NS mergers. 
Very recent calculations [70l. [7ll [72j indicate radiative 
energies ranging from ~ 10~ 4 to ~ 10~ 2 of the total 
mass, where the range is likely to reflect the very differ- 
ent initial conditions assumed in the simulations. While 
there are no short-duration GRBs in the S4 sample, we 
can use typical upper limits on h Taa from Table I VII as 
an indication of sensitivity. For example a 1.4M Q NS 
plus 10M Q BH binary system would have merger GW 
emission at frequencies starting at about 400 Hz. If this 
system were to radiate 1% of its rest energy into gravita- 
tional radiation at 400 Hz, the distance sensitivity would 
be D ~ 5 Mpc. The search would also be sensitive to the 
inspiral emission from this system at lower (~ 200 Hz) 
frequency. 



D. Prospects 

Here we discuss the future prospects for science run 
S5 and beyond. At the sensitivity for science run S4, 
the prospects for detection are clearly dominated by the 
possibility of a nearby GRB. While this distance scale is 
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guided by the discussion above, we are prepared to be 
surprised by new mechanisms for GW emission. Never- 
theless, we expect detection of individual GRBs to de- 
pend in no small part on the appearance of a "special" 
event. Thus, a data sample which includes a large num- 
ber of GRBs is especially important. For science run S5, 
the LIGO detectors will be operating at design sensitiv- 
ity and fully coincident with Swift operation. This should 
yield over 100 GRBs, including some with redshift mea- 
surements. And clearly, the search radius will increase in 
proportion to improvements in the LIGO strain sensitiv- 
ity. 

The results pertaining to the GRB population obtained 
in Section IVIII will certainly improve for the S5 run and 
in future observations with Advanced LIGO. To make 
an estimate, we look at the various factors involved in 
Eq. [201 for the upper limit on £?gw- As a reference, we 
use the limit obtained here using all S2, S3 and S4 GRBs. 
Since most factors in Eq. [20] come as squares, moderate 
improvements in each has a significant overall effect. 

Since the direction to each GRB will be known, it may 
be possible to select a subsample of, say, 35 GRBs from 
the sample in S5 (i.e., about the same number as the 
whole of S2, S3 and S4) such that (F avc ) ~ 1/V2, the 
maximum possible. Further, assume that we use Hl-Ll 
crosscorrelations. Figure [TU] shows the confidence belt 
for the case of 35 optimally located GRBs and a pair of 
identical detectors. One can expect to get an upper limit 
of ~ 10 on po with this curve, which is a factor of ~ 3.5 
better than the current limit on pq. 

Without altering other parameters of the analysis, 
therefore, we can expect 3.5 2 or, in round numbers, a 
factor of ~ 10 improvement in the upper limit on £gw 
for S5. Additional improvements are possible by impos- 
ing a cut based on measured redshifts, in addition to the 
cut on sky positions, and by reducing the search interval 
from the current value of 180 seconds. Looking beyond 
S5, the most obvious source of improvement would be 
the > 10 factor of improvement in the strain noise level 
when Advanced LIGO comes online around the middle 
of the next decade. This translates into an additional 
factor of > 100 reduction in the upper limit. When Ad- 
vanced LIGO comes online, there may be a worldwide 
network of GW detectors of comparable sensitivity. Be- 
sides allowing a more uniform sky coverage, resulting in 
a larger sample of GRBs with optimal orientation, net- 
work analysis methods [73|, [74], [75| that make more opti- 
mal use of data from multiple detectors can be used to 
increase the base sensitivity of the method. Finally, with 
enough GRBs, we could separately analyze the class of 
long and short duration bursts. Since the most proba- 
ble redshift for short-duration GRBs is expected to be 
inherently smaller, we could obtain significantly tighter 
constraints on the energy emitted in gravitational waves 
from this class of GRBs. 

The discussion above was confined to a particular 
model for GRB redshift distribution and GW emission. 
Further work is needed to develop more general analysis 



methods that can be applied to a wider variety of models 
and that take better account of prior information from 
existing observations. 



IX. SUMMARY AND CONCLUSION 

We searched for gravitational-wave bursts, target- 
ting short GW signals with durations from ~ 1 ms to 
~ 100 ms, associated with 39 GRBs that were detected 
by gamma-ray satellite experiments while the S2, S3, and 
S4 science runs of the LIGO experiment were in progress. 
To take into account the unknown onset time of the GW 
signal relative to the GRB trigger time, the search cov- 
ered 180 seconds of data surrounding the GRB trigger 
times. These 180-second data segments from the dif- 
ferent IFOs were crosscorrelated to probe for correlated 
signals. We searched for an association on an individual- 
GRB basis, and also applied different statistical tests to 
search for the cumulative effect of weak GW signals. We 
found no evidence for gravitational-wave burst emission 
associated with the GRB sample examined using the dif- 
ferent search methods. 

Using simulated Q = 8.9 sine-gaussian waveforms and 
the direction-dependent antenna response of the inter- 
ferometers to a GW source, we obtained upper limits on 
the root-sum-square amplitude of linearly polarized and 
circularly polarized gravitational waves from each of 22 
GRBs with well-localized positions. Associating these 
limits with the energy radiated by the GRB sources into 
gravitational radiation is inherently speculative at this 
stage of development of the field and depends crucially 
on the astrophysical scenario one adopts for the GRB 
progenitors. The most favorable cases considered here 
suggest that the LIGO sensitivity for run S4 would allow 
sensitivity to a solar mass-equivalent of radiated GW en- 
ergy to distances of tens of Mpc. 

The sample of GRBs was combined to set an upper 
limit on the GW energy emitted using a simple standard 
candle model and a theoretical redshift distribution of 
GRBs. Although the upper limit obtained is not astro- 
physically important, a straightforward and realistic ex- 
trapolation to future observations suggests that this limit 
can be improved by orders of magnitude. It may be pos- 
sible to set a sub-solar mass limit when Advanced LIGO 
comes online. This would put us in an astrophysically in- 
teresting regime since at least one model [761 ] predicts an 
energy loss of 0.2 solar masses for long-duration GRBs. 

It is opportune that Swift will be operating and de- 
tecting GRBs at the time when the fifth science run of 
LIGO, S5, will be in progress. The goal for the S5 run 
is to collect one year of coincident LHO-LLO data at the 
design sensitivity. Given the Swift GRB detection rate, 
we anticipate an S5 sample of more than 100 GRB trig- 
gers that can be used to further probe for gravitational 
radiation associated with GRBs. It is hoped that a large 
GRB sample will increase the chances for finding such an 
association. 
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TABLE III: S4 90% upper limits on hrss of Q = 8.9 linearly polarized sine-gaussians, in units of 10 21 Hz 1//2 ; 25-ms crosscorrelation length. 
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TABLE IV: S3 90% upper limits on hrss of Q = 8.9 linearly polarized sine-gaussians, in units of 10~ 20 Hz~ 1/2 ; 25 -ms crosscorrelation length. 
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TABLE V: S2 90% upper limits on hrss of Q = 8.9 linearly polarized sine-gaussians, in units of 10~ 19 Hz~ 1/2 ; 25 -ms crosscorrelation length. 
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7.1 


2.5 


1.1 


2.2 


0.69 


0.58 


0.71 


1.0 


1.1 


1.3 


1.6 


2.9 


3.1 


3.8 


6.0 


5.6 


030323a 


5.1 


3.1 


6.4 


2.5 


1.7 


2.9 


1.1 


0.99 


1.5 


1.7 


2.3 


3.3 


2.6 


6.1 


7.2 


6.0 


11.4 


13.4 


030323b 


4.6 


1.8 


5.2 


1.7 


0.94 


1.8 


0.64 


0.45 


0.81 


0.92 


0.82 


1.5 


1.3 


1.8 


2.4 


3.0 


3.5 


4.8 


030324 


9.2 






4.7 






1.6 






2.0 






3.3 






7.9 






030325 


2.8 


1.7 


3.0 


1.3 


0.80 


1.5 


0.55 


0.48 


0.76 


0.89 


1.0 


1.5 


1.3 


2.0 


2.4 


3.2 


4.9 


5.3 


030326 


10.2 


3.9 


9.6 


4.4 


2.1 


3.7 


1.4 


0.94 


1.2 


2.0 


1.6 


1.9 


3.1 


3.4 


3.1 


8.4 


8.1 


6.3 


030329a 


4.6 






2.4 






1.1 






1.8 






3.0 






7.6 






030329b 


2.8 






1.1 






0.31 






0.55 






0.89 






2.0 






030331 




3.4 






1.6 






0.85 






2.0 






3.4 






8.0 




030405 


2.1 


1.4 


3.1 


1.0 


0.80 


1.3 


0.34 


0.42 


0.51 


0.59 


0.76 


0.97 


0.87 


2.0 


2.2 


2.0 


4.8 


4.5 


030406 




1.2 






0.67 






0.42 






0.77 






1.7 






4.4 




030413 






1.6 






0.85 






0.50 






0.89 






2.3 






4.4 


030414 


1.4 






0.91 






0.32 






0.39 






0.70 






1.6 







TABLE VI: S4 90% upper limits on hrss of Q = 8.9 circularly polarized sine-gaussians, in units of 10~ 21 Hz~ v ; 25 -ms crosscorrelation length. 





100 Hz 


150 Hz 


250 Hz 


554 Hz 


1000 Hz 


1850 Hz 


GRB 






































date 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


050223 


1.6 






1.1 






1.2 






2.0 






3.5 






6.7 






050306 


2.2 


1.4 


2.6 


1.5 


1.1 


1.8 


1.6 


1.4 


2.0 


2.6 


2.6 


3.3 


4.5 


5.0 


6.2 


8.5 


14.2 


17.6 


050318 


2.2 


2.2 


3.1 


1.6 


1.5 


2.2 


1.6 


1.9 


2.4 


2.6 


3.5 


4.0 


4.6 


6.1 


6.9 


8.8 


11.1 


13.1 


050319 


1.8 


1.8 


2.3 


1.4 


1.3 


1.6 


1.5 


1.7 


1.8 


2.4 


3.1 


3.0 


4.3 


5.5 


5.2 


8.2 


10.0 


9.9 



TABLE VII: S3 90% upper limits on hrss of Q = 8.9 circularly polarized sine-gaussians, in units of 10 1 Hz ' ; 25-ms crosscorrelation length. 



GRB 
date 


100 Hz 


150 Hz 


250 Hz 


554 Hz 


1000 Hz 


1850 Hz 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


031108 
031109a 
031220 


19.0 
14.7 
14.4 






11.3 
8.8 
10.1 






10.9 
8.5 
8.9 






12.5 
10.6 
10.8 






20.4 
17.3 
18.4 






51.5 
42.2 
42.7 







TABLE VIII: S2 90% upper limits on hrss of Q = 8.9 circularly polarized sine-gaussians, in units of 10~ 20 Hz _1/2 ; 25 -ms crosscorrelation length. 





100 Hz 


150 Hz 


250 Hz 


554 Hz 


1000 Hz 


1850 Hz 


GRB 






































date 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


030226 


22.2 


11.0 


18.0 


9.2 


5.0 


6.9 


2.9 


2.1 


1.9 


3.7 


3.3 


2.6 


7.1 


7.1 


4.1 


20.3 


20.3 


7.2 


030320a 


21.9 


7.0 


26.6 


7.3 


3.6 


7.9 


2.0 


1.9 


2.2 


2.9 


3.3 


4.1 


4.6 


9.5 


10.1 


10.7 


17.3 


16.1 


030323a 


16.1 






7.9 






3.6 






5.6 






7.9 






18.5 






030323b 


13.4 


4.9 


15.5 


4.9 


2.5 


5.1 


1.8 


1.2 


2.3 


2.7 


2.4 


4.0 


3.7 


5.1 


6.6 


8.5 


9.2 


12.3 


030324 


28.0 






13.3 






4.3 






5.6 






9.4 






22.2 






030325 


9.0 


4.3 


9.5 


4.0 


2.0 


4.2 


2.0 


1.2 


2.4 


3.1 


2.8 


4.4 


4.3 


5.3 


6.7 


10.2 


12.2 


15.0 


030326 


29.7 


15.1 


39.9 


12.4 


8.1 


14.9 


4.0 


3.5 


4.8 


5.8 


5.8 


7.6 


9.6 


12.1 


11.7 


24.2 


25.8 


19.7 


030329a 


13.8 






7.3 






3.3 






5.1 






8.2 






21.6 






030329b 


8.8 






3.2 






0.90 






1.5 






2.4 






5.9 






030331 




7.1 






3.5 






1.8 






4.3 






7.3 






17.4 




030405 


6.2 


3.4 


8.2 


2.9 


2.0 


3.4 


0.99 


1.1 


1.3 


1.6 


2.0 


2.5 


2.5 


5.1 


5.4 


5.9 


11.3 


10.7 


030406 




2.8 






1.5 






0.90 






1.8 






4.0 






10.0 




030413 






4.1 






2.2 






1.3 






2.4 






6.0 






11.0 


030414 


4.1 






2.6 






0.82 






1.1 






1.9 






4.6 







TABLE IX: S4 90% upper limits on hrss of Q = 8.9 linearly polarized sine-gaussians, in units of 10 21 Hz 1 ' 2 ; 100- ms crosscorrelation length. 





100 Hz 


150 Hz 


250 Hz 


554 Hz 


1000 Hz 


1850 Hz 


GRB 






































date 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


050223 


5.6 






4.1 






4.8 






8.0 






14.5 






30.9 






050306 


6.9 


6.7 


12.6 


4.9 


5.8 


9.1 


5.6 


7.6 


10.4 


9.1 


13.8 


17.3 


16.0 


28.0 


34.0 


30.0 


74.1 


91.8 


050318 


7.4 


9.7 


12.5 


5.9 


7.4 


10.3 


6.4 


9.9 


11.8 


10.7 


17.5 


17.9 


18.4 


33.2 


34.1 


33.3 


63.4 


64.5 


050319 


5.5 


6.0 


9.6 


4.6 


4.6 


7.2 


5.2 


6.5 


8.4 


8.8 


11.4 


14.4 


15.2 


21.3 


25.1 


30.1 


34.7 


48.3 



TABLE X: S3 90% upper limits on hrss of Q = 8.9 linearly polarized sine-gaussians, in units of 10~ 20 Hz~ 1/2 ; 100 -ms crosscorrelation length. 





100 Hz 


150 Hz 


250 Hz 


554 Hz 


1000 Hz 


1850 Hz 


GRB 






































date 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


031108 


6.0 






3.6 






3.8 






4.5 






7.9 






20.1 






031109a 


4.4 






2.7 






2.9 






3.5 






6.1 






15.1 






031220 


5.0 






3.0 






3.0 






4.1 






7.0 






15.8 







TABLE XI: S2 90% upper limits on hrss of Q = 8.9 linearly polarized sine-gaussians, in units of 10~ 19 Hz~ 1/2 ; 100 -ms crosscorrelation length. 





100 Hz 


150 Hz 


250 Hz 


554 Hz 


1000 Hz 


1850 Hz 


GRB 






































date 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


030217 






4.0 






2.0 






0.94 






1.5 






4.1 






9.5 


030226 


7.3 


3.1 


5.3 


3.2 


1.5 


2.1 


1.1 


0.65 


0.62 


1.4 


1.0 


0.85 


2.6 


2.4 


1.4 


7.1 


6.5 


2.7 


030320a 


6.7 


2.3 


6.8 


2.5 


1.3 


2.3 


0.76 


0.67 


0.70 


1.0 


1.2 


1.4 


1.6 


3.6 


3.5 


4.1 


7.2 


5.8 


030323a 


5.3 


2.7 


5.6 


3.0 


1.5 


2.5 


1.2 


0.86 


1.4 


1.8 


2.2 


3.0 


2.7 


5.5 


7.0 


6.4 


10.0 


12.4 


030323b 


5.1 


1.8 


4.9 


2.0 


0.95 


1.7 


0.77 


0.47 


0.79 


1.1 


0.90 


1.6 


1.6 


1.9 


2.5 


3.9 


3.7 


5.0 


030324 


8.7 






4.6 






1.5 






2.0 






3.7 






8.0 






030325 


2.9 


1.5 


3.4 


1.4 


0.78 


1.6 


0.63 


0.46 


0.90 


1.0 


1.00 


1.9 


1.5 


1.9 


2.9 


3.7 


4.6 


6.6 


030326 


9.0 


3.0 


7.4 


4.2 


1.8 


3.1 


1.3 


0.81 


0.98 


1.9 


1.5 


1.8 


3.7 


3.1 


2.9 


8.6 


6.8 


5.7 


030329a 


4.4 






2.5 






1.2 






2.1 






3.0 






8.6 






030329b 


2.6 






1.2 






0.34 






0.56 






0.94 






2.2 






030331 




3.5 






1.7 






0.97 






2.1 






4.1 






10.3 




030405 


2.3 


1.2 


2.6 


1.3 


0.76 


1.1 


0.46 


0.40 


0.47 


0.73 


0.73 


0.90 


1.2 


1.8 


1.9 


2.7 


4.4 


4.0 


030406 




1.2 






0.73 






0.45 






0.87 






1.9 






5.0 




030413 






1.7 






0.94 






0.61 






1.1 






2.9 






5.4 


030414 


1.3 






0.89 






0.30 






0.43 






0.74 






1.7 







TABLE XII: S4 90% upper limits on hrss of Q = 8.9 circularly polarized sine-gaussians, in units of 10~ 21 Bfz 7 ; 100 -ms crosscorrelation length. 





100 Hz 


150 Hz 


250 Hz 


554 Hz 


1000 Hz 


1850 Hz 


GRB 






































date 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


050223 


1.7 






1.3 






1.5 






2.4 






4.4 






8.3 






050306 


2.0 


1.5 


2.6 


1.5 


1.2 


1.9 


1.7 


1.7 


2.2 


2.8 


3.1 


3.7 


4.9 


6.0 


7.0 


9.3 


16.3 


19.1 


050318 


2.2 


2.1 


2.8 


1.7 


1.6 


2.2 


1.9 


2.1 


2.4 


3.0 


4.0 


4.2 


5.5 


6.9 


7.4 


10.3 


12.7 


14.0 


050319 


1.7 


1.6 


2.5 


1.4 


1.3 


1.9 


1.6 


1.7 


2.2 


2.6 


3.2 


3.8 


4.7 


5.7 


6.7 


9.1 


10.3 


12.8 



TABLE XIII: S3 90% upper limits on hrss of Q = 8.9 circularly polarized sine-gaussians, in units of 10~ 21 Hz~ 1/2 ; 100 -ms crosscorrelation length. 



GRB 
date 


100 Hz 


150 Hz 


250 Hz 


554 Hz 


1000 Hz 


1850 Hz 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


031108 
031109a 
031220 


18.4 
13.5 
12.1 






11.5 

8.5 
9.4 






11.8 
8.7 
8.8 






14.0 
11.3 
11.6 






23.2 
19.0 
20.5 






61.0 
47.6 
49.1 







TABLE XIV: S2 90% upper limits on hrss of Q = 8.9 circularly polarized sine-gaussians, in units of 10~ 20 Hz^ 1/2 ; 100 -ms crosscorrelation length. 





100 Hz 


150 Hz 


250 Hz 


554 Hz 


1000 Hz 


1850 Hz 


GRB 






































date 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


H1-H2 


Hl-Ll 


H2-L1 


030226 


22.1 


9.6 


16.7 


9.8 


4.7 


6.1 


3.1 


2.0 


1.9 


4.1 


3.2 


2.6 


7.8 


7.2 


4.5 


21.2 


19.6 


8.0 


030320a 


21.0 


7.4 


24.7 


7.7 


4.1 


7.6 


2.2 


2.2 


2.4 


3.3 


4.0 


4.5 


5.3 


11.4 


10.8 


12.3 


21.7 


18.3 


030323a 


16.7 






8.8 






4.0 






6.3 






9.3 






21.8 






030323b 


14.8 


4.8 


14.4 


5.9 


2.6 


4.8 


2.2 


1.3 


2.4 


3.3 


2.5 


4.4 


4.7 


5.3 


7.4 


10.9 


10.4 


14.1 


030324 


27.0 






13.9 






4.7 






6.3 






10.7 






24.7 






030325 


9.7 


3.7 


9.9 


4.6 


2.0 


4.5 


2.2 


1.2 


2.8 


3.5 


2.6 


5.4 


5.2 


5.1 


8.4 


12.3 


12.5 


19.1 


030326 


28.3 


11.0 


28.6 


13.0 


6.3 


10.9 


4.3 


2.9 


3.7 


6.3 


5.0 


6.2 


10.6 


10.6 


10.2 


26.4 


23.2 


18.3 


030329a 


13.7 






7.8 






3.6 






5.8 






9.5 






24.8 






030329b 


8.1 






3.3 






1.0 






1.7 






2.8 






6.7 






030331 




7.4 






3.7 






2.1 






4.9 






8.6 






20.6 




030405 


7.1 


3.1 


6.8 


3.7 


1.9 


2.9 


1.3 


1.1 


1.2 


2.1 


1.9 


2.3 


3.3 


5.0 


5.2 


7.8 


11.6 


10.4 


030406 




2.8 






1.7 






1.0 






2.0 






4.6 






11.6 




030413 






4.3 






2.4 






1.5 






2.8 






7.3 






13.5 


030414 


4.1 






2.7 






0.91 






1.3 






2.2 






5.2 







